Spin torque antiferromagnetic nanooscillator in the presence of the magnetic noise 
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Spin-torque effects in antiferromagnetic (AFM) materials are of great interest because of the 
possible applications as high-speed spintronic devices. In the present paper we analyze the statistical 
properties of the current-driven AFM nanooscillator that result from the white Gaussian noise of 
magnetic nature. Starting from the peculiarities of deterministic dynamics we derive Langevin and 
Fokker-Planck equations in energy representation for two normal modes. We find the stationary 
distribution function in the subcritical and overcritical regimes and calculate the current dependence 
of average energy, energy fluctuation and their ratio (quality factor). It is found that one of the 
modes (noncritical) shows the Boltzmann statistics with current-dependent effective temperature in 
all range of current values. The effective temperature of the other (soft) mode critically depends 
on current in the subcritical region. Distribution function of soft mode follows the Gaussian law 
above the generation threshold. In the overcritical regime the full average energy and quality factor 
grow with current value that makes AFM nanooscillators promising candidates as active spintronic 
components. 
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I. INTRODUCTION 



Nowadays the spin-polarized current is widely used for manipulation of nano-magnetic structures. Corresponding 
physical mechanism is based on the so-called spin-transfer-torque (STT) effect, which was predicted by Slonczewski 
and Berger [1, 2]: a spin-polarized current may transfer angular momentum to a free ferromagnetic (FM) layer 
resulting in a macroscopic torque on the latters magnetization. In small FM particles (nanopillars or nanocontacts) 
in which the magnetization can be assumed to be spatially uniform, the STT results in the steady rotation of the 
magnetization (see, e.g. [3-6]). This effect is interesting for applications: as a tempting alternative to the electronic 
(us spintronic) diode, as radio-frequency devices used for telecommunications, as timiming mechanisms etc. 

Recently it was shown [7] (see also [8]) that the STT effect should take place also in antiferromagnetic (AFM) 
materials and, in analogy with FMs, should give rise to steady rotation of the Neel (or AFM) vector. Current- 
controlled AFM nanoparticles are promising candidates for spintronic devices because of high working frequencies that 
fall into 0.1-1 THz range (for comparison, typical frequencies of FM nanooscillators are 1-50 GHz [9, 10]). However, 
the central problem for practical applications of nanooscillators is the reduction and control of their linewidth and 
as a consequence, improvement of their quality factor. Thus, understanding of the stochastic processes (such as, e.g. 
thermal noise) that set condition for linewidth is increasingly important for operation of STT devices. 

Thermal noise in FM particles was a matter of study starting from the pioneering work by W. F. Brown [11], 
where the noise has been described as a stochastic magnetic field acting on the magnetization and corresponding 
Fokker-Planck equation was derived. Apalkov and Visscher [ 2] generalized the nonlinear Fokker-Planck equation 
to the case of Slonczewski STT. Since then noise properties of the FM-based nonlinear oscillator in the presence of 
spin-polarized current were studied both experimentally and theoretically [13-16]. 

Investigation of noise in AFM systems is much more complicated problem compared to that in FMs due to i) the 
larger number degrees of freedom; ii) Newtonian-like (vs precession-like in FMs) dynamics of the Neel vector. In some 
cases, e.g. the description of superparamagnetism or switching processes in AFM nanoparticles, the problem can be, 
in principle, reduced to the description of fluctuations of weak FM moment which inevitably arises due to the surface 
effects or Dzyaloshinskii-Moriya interactions [17-20]. 

However, those peculiarities of AFM dynamics that are related with the exchange coupling between the magnetic 
sublattices, like magnetoelastic effects [21], spin-wave spectra [22], STT phenomena [7] etc, should be described with 
due account of the variables inherent to AFM ordering. Corresponding approach based on the dynamics of the Neel 
vector was recently developed [23] for the case of a collinear compensated AFM with two magnetic sublattices. In 
the present paper we apply this approach to the analysis of the efficiency of AFM nanooscillator conducted by spin- 
polarized current. We assume that the thermal noise in AFM particle arises from fluctuations of the random magnetic 
field. Following the method of slow and fast variables [24] and energy representation for non-equilibrium steady state 
[25] we formulate the Fokker-Planck equation for energy distribution and study the linewidth of spin-torque AFM 
nanooscillator depending on temperature and current. 



II. MAGNETIC DYNAMICS OF ANTIFERROMAGNETIC PARTICLE IN THE PRESENCE OF 

SPIN-POLARIZED CURRENT 

The system under consideration (see figure 1 a) is an AFM nanoparticle with the elliptic crossection. The size 
of the particle is large enough to ensure the AFM ordering and is small enough (below the domain wall thickness, 
so-called macrospin approximation) to neglect the space variation of the magnetic properties. Hard FM (polarizer) 
and thin nonmagnetic layers shown in figure la) are the ancillary elements that deliver the spin-polarized current to 
AFM. We assume that the temperature T of the system is kept constant thus neglecting the thermal (Joule) heating 
of the system. 

For the sake of simplicity we consider a collinear AFM with two equivalent magnetic sublattices Mi and M2 
and disregard weak FM moment that can arise either from the intrinsic properties of material or due to the surface 
effects. The coupling between the sublattices (characterized with the spin- flip field He) is assumed to be much 
stronger than other fields (including anisotropy and all the external fields). To this end, the magnetic state of such 
an AFM nanoparticle is unanimously described by the only one, AFM vector, L = Mi — M2 of the fixed length 
|L| = 2Mq. Corresponding dynamic equations for L can be obtained within the standard Lagrange technique, the 
Lagrange function being [26] 



Cafm = ^L 2 + 7 m L [L • (L x H)] - ™ an (L) + ^^(L x H) 2 , (1) 

where H is an external magnetic field and w aa (L) is the energy of magnetic anisotropy that forms a potential well 
for AFM vector, 7 is the gyromagnetic ratio. The value tul = 1 / (2^ 2 MqH e) plays a role of the "inertia mass" of the 
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Figure 1. Dynamics of AFM vector in the presence of spin-polarized current, a) AFM particle of the elliptic cross-section is 
placed between two electrodes, the top electrode being FM (polarizer). The spacer between FM and AFM is nonmagnetic to 
avoid direct exchange coupling between the magnetic layers, b) Typical trajectories of L vector for circularly polarized modes. 
Different areas correspond to different equilibrium orientations (4 of six possible are shown), c) Projection of trajectories in L 
space to the xy plane. Three inner circles correspond to the normal modes with different amplitude in the vicinity of equilibrium 
state L || 2;. The last but one circle corresponds to separatrix and the outer circle is a steady state trajectory, d) Trajectory in 
the overcritical regime (J = 1.2J cr ). In the initial state L||p curr ||z. 

Newtonian "material point" with "radius-vector" L. 

For the definiteness, in what follows we assume that the AFM layer has a cubic magnetic anisotropy modeled with 
the following expression 



8M$ 



(Li + Lt + Li), 



(2) 



where the orthogonal axes x, y and z coincide with the easy directions for the Neel vector, H an -C He is anisotropy 
fields. 

Dissipation is modeled with the Rayleigh function which in the presence of spin-polarized current J takes a form 



R-afm = 7afm"^lL 



crJ 

2jM 



(L x L) 



(3) 



Here the first term models the internal damping, damping coefficient 27afm is, in fact, the AFMR linewidth, the 
constant a = Ti^e/ [leM^u afm) is proportional to the efficiency e of the spin transfer processes, vafm is the volume 
of AFM nanoparticle, fi is the Plank constant, e is the electron charge. Unit vector p cuxr is parallel to the direction 
of the current spin polarization. 

Deterministic behaviour of such a system in the presence of current J polarized along one of the easy axes, 
Pcurr||z, was analyzed in details in reference [7]. Some characteristic features important for further consideration 
are summarized in figure l(b-d). In the so-called subcritical regime, |J| < J cr j t = 27afm^afmr/(7C-H£;) (where 
Oafmr = iVHeHIh is AFMR frequency) the AFM vector L has three equilibrium orientations close to x, y and z 
directions that define six (corresponding to L and — L) basins of finite motion in the phase space. Small deflection 
from the easy axis in equilibrium is due to the current-induced torque. Typical phase trajectories in the vicinity of 
equilibrium points in nondissipative approximation correspond to the clockwise/counterclockwise rotations of AFM 
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vector (circles in figure lb) and could be associated with the circular polarized normal modes with the frequencies 
±^afmr- Spin-polarized current acts as a negative damping for one of the modes ("soft" mode) and as a positive 
damping for the other. However, below the critical current the negative damping is suppressed by the internal losses 
and the real phase trajectories are the twisted spirals. 

In the overcritical regime (|J| > J cr it) the stable (non-equilibrium) state corresponds to steady rotation of AFM 
vector in xy plane (perpendicular to p CU rr) with the current-dependent frequency lu — ( J/ J cr it)^AFMR (the outer 
circle in figure lc). The direction of the rotation coincides with the direction of the soft mode. In the steady state 
the internal losses are exactly compensated with the current-induced energy pumping. This regime can be associated 
with power generation. 

Trajectory of AFM vector presented in figure Id) shows the typical current-induced transition from the equilibrium 
state L||z to the steady precession in xy plane. The prominent feature of this motion is the presence of two well 
separated time scales: fast rotation around z axis with the frequency cx ^afmr and slow, with the characteristic time 
oc 1/7afm variation of the polar angle from (L||z) to tt/2 (L _L z). Deviation from this scenario takes place only 
in the close vicinity of the separatrix (see figure lc, next to last circle) where the rotation frequency substantially 
diminishes. 

Thus, in the vicinity of equilibrium and stationary steady states the current-induced dynamics of AFM is charac- 
terized with the set of slow and fast variables that allows substantial simplification of the description. 



III. LANGEVIN DYNAMICS AND FOKKER-PLANCK EQUATIONS IN ENERGY REPRESENTATION 

While the state of FM nanoparticle is described with only two variables in configuration space, the dynamics of 
AFM nanoparticle needs for description at least four (with account of normalization condition |Mi| = IM2I = Mo 
far below the Neel temperature) variables in the phase space, namely, generalized coordinates L and corresponding 
generalized momenta Pl. This results in a rather complicated Langevin equations in the phase space [23]: 

L = P L /m L - 7 L x h, (4) 

Pl = F L + F diss - 7 (P L - 27 AFM m i L) x h, 

where = — (9w an (L)/<9L is the potential (gradient) force, and the dissipative force Fdi ss is given by the following 
expression 



, _ <97?-AFM 

diss = 



<9L 



= -27afmPl - Pcurr x L. (5) 



The thermal noise in equation (4) is modeled with the random magnetic field h(t) which we assume to be the white 
Gaussian noise with 

(h(t)) = 0, {h^hkih)) = 2D6 jk 6(h - t 2 ), (6) 

where 2D represents the intensity of thermal fluctuations. 

However, the number of the effective phase variables in Langevin, (4), and corresponding Fokker-Planck equations 
can be reduced if we take into account the above mentioned peculiarities of the dynamics in the presence of spin- 
polarized current. 

First, in the vicinity of equilibrium and stationary states the motion of AFM vector is finite and can be decomposed 
to a linear combination of two independent (normal) modes (see figure 2). This means that one can use the set 
of the canonically conjugated variables action-angle, I±,f± (± correspond to the clock/counterclockwise rotations 
in configuration space) instead of coordinates and momenta. Second, in a low-temperature approximation, when 
the temperature T is much less than the energy barrier between the different equilibrium/stationary states, all the 
essential phase trajectories for each mode have the same rotation/oscillation frequency, lo±. Thus, instead of action 
one can use the energy of the mode, E± = I±u>±, as a canonical variable. At last, two time scales allow one to exclude 
the angle variables ip± by averaging over the period of rotations. 

As a result, Langevin equations (4) could be rewritten in terms of the averaged energies E± as follows: 



dE± 
dt 



' diss + 7 



2 7 afmPl - -TjjT- ) ■ L x h 



(7) 



where the overline means averaging over the period of rotation and the summands in r.h.s. of equation (7) should be 
expressed in terms of E± . Here the averaged energy is 

/>2?r / LJ± 

E ± = ^ [P|/(2m L ) + Wan (L)] dt (8) 




Figure 2. Normal circularly polarized modes in the presence of spin-polarized current in a) subcritical and b) overcritical 
regimes. Typical trajectories lay over the sphere |L| = 2Mq in configuration space. 



The explicit closed form of equation (8) will be obtained in next subsections for the limiting cases of subcritical and 
overcritical regimes. 



A. Subcritical regime, J < J cr 

Let us consider the basin of states in the vicinity of equilibrium point L||p curr ||z and chose L x ,L y <C 2Mq as 
generalized coordinates, L z ps 2M . For the normal modes the time dependence of the dynamic variables is given by 
the expressions 

L x + iL y = 2M c±e iu±t , P Lx + iP Ly = 2iM m L uj±c±e lul±t , lu± = ±n AFMR . (9) 

The averaged energy is related with the amplitude c± as follows: E± — 2E c±, where the value E = 2MqQ\ fmr itil = 
Mo-ff a n defines the characteristic energy scale for AFM nanoparticle. 

Substituting expressions (9) into (5) and (7) we arrive to the system of two independent Langevin equations: 

= -27afm (l ± -y— ) £± + 2j\/£± (±h x sin fi A FMRi - h y cos Qafmr*) (10) 
at \ J cr it / 

+ 7 ^ 7AFM cos Oafmr* ± /lySinOAFMR*) T£±h z ] , 

s 'afmr. 

where e± = E±/Eq is a dimensionless energy. 

As it is seen from equation (10), the clock/counterclockwise modes (figure 2a) interact with the current in different 
ways. If J > 0, the effective damping of the first mode (with the energy E+) increases and that of the second (with 
the energy _E_) decreases, due to the action of spin-polarized current. 

Equations (10) are the typical Langevin equations in energy representation considered in details in [25]. The first 
summand in the r.h.s. describes the rate of the direct energy exchange which depends upon the current value J. All 
but first summands in the r.h.s. of (10) account for system-environment interaction induced by field fluctuations. 
The diffusion functions (coefficients before hj) depend on energy and thus correspond to multiplicative noise. Last 
two terms (in square brackets) are multiplied by small factor 7afm / ^afmr <C 1 and will be neglected for the sake of 
simplicity. 

In the accepted approximation of nonintcracting modes the distribution function in phase space, /(L, Pj,;£), can 
be factorized as follows: /(L,Pl;£) =>■ f+(E+;i)f-(E-;t). The Fokker-Planck equations for f± are then deduced in 
a standard manner from (10) in Stratanovich convention and take a form: 



df±(E±) _ 8 



dt dE ± 



1 2 DE ^E~Z 7 ^-^eZ+2 1afm (l±^-)Ei 

OE± V J crit 



f±(E±)\. (11) 



From the stationary solution of (11) one gets the AFM probability distribution function f(E+,EJ) — f + (E + )f-(E-)y/E~ 
(with due account of Jacobian): 



n r r , r I 2 7 A FM 

f(E+,E-) = fo exp <^ — — 



(12) 
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where fo is a normalization constant and the diffusion coefficient is related with the temperature T through the 
fluctuation-dissipation theorem as follows: 

D = *[™LT. (13) 
T e q 



B. Overcritical regime, J > J cr 

In the overcritical regime the motion can be decomposed (like it was done in the reference [27] for FM) into steady 
rotation of AFM vector in xy plane with the frequency w_ and small meridian oscillations of L vector with the 
frequency lo + (see figure 2b). 

The "+"-mode is parametrized as follows 

L z = ce iu +\ P Lz = 2iM m L uj + ce l "+ t , 



3 J 2 /3 J 2 \ 

u+ = r»AFMRW - + -p, E+ = 2E (- + - c 2 . (14) 

Corresponding Langevin and Fokker-Planck equations are analogous to (10) and (11). 

Parametrization of "—"-mode coincides with that in equation (9) with c_ = 1. In this particular case the proper 
dynamic variable is action Z_ = ^irM^rriLUj^. Corresponding Langevin equation takes a form 

^- = -2 7 Afm(/- - 87rAf 2 m LW l 0) ) - jl_h z , (15) 

where wf ' = — (J/ J C r)^AFMR is the deterministic frequency of steady rotation. Transition to energy representation 
can be obtained directly from (15) by substitution J_ = AttMq J 2m^E^ : 

dE- {„ J 



dt \ J c , 

Corresponding Fokker-Planck equation is then deduced as following: 



4 7AF M E- - — VE E- - 2 1 E^h z . (16) 



df-{E-) d 



<kj 2 DE_-^E_ + 4 7AFM \E- - j-^/E E_ 



dt dE- 

Ultimately, the stationary AFM distribution function in the overcritial regime takes a form: 



/-(£-) • (17) 



f(E+,E.) = / exp 



4£+ (E--E J 2 /J 2 r ) 2 



[3 + 4(J/J crit ) 2 ]T 2TE Q J 2 /J 2 r 
where we have taken into account the relation (13) and assumed that T <C Eq. 



(18) 



IV. DISCUSSION 



In the previous section we derived the stochastic equations for the current-controlled AFM nanoparticle in the 
low temperature approximation. The energy scale of the system is given by the magnetic anisotropy energy Eq, so, 
the validity of equations obtained is limited by inequality T ^ Eq. Taking into account that characteristic value of 
magnetic anisotropy for the AFMs with high Neel temperature is 10 3 4- 10 4 J/m 3 (see, e.g. [28] for Mns2Niig) and 
typical size of nanoparticle is 50 x 50 x 5 nm 3 one gets the following estimation: Eq oc 10~ 20 -j- 10 -19 J. Thus, the 
proposed model can be applied up to the room temperature, Trt = 4 • 10 -21 J. 

We found the stationary solutions of Fokker-Planck equations that describe the distribution function in the vicinity 
of equilibrium, (12), and nonequilibrium steady, (18), states where the swapping between these two states can be 
neglected. 

For the noncritical ("+") mode the probability function ,f(E + ) (see (12) and (18)) follows the Boltzmann law with 
the current-dependent effective temperature 

T->(+) — T 1 / (1 + J/Jcrit) 1 J < >^crit: / 1Q \ 

- 1 \ 0.75 + J 2 /J 2 rit , J> J crit . • (1Jj 
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Figure 3. Energy distribution function for different current values J/J CT (shown with numbers near the curves) in the subcritical 
(J < J cr ) and overcritical (J > J cr ) regimes. Energy is measured in the dimensionless units, E/T. Vertical dashed line shows 
the energy Eq — 10T in the overcritical regime. 



The "soft" ("— ") mode shows the Boltzmann-like distribution in the subcritical regime (J < J cr ;t) which changes 
to the Gauissian-like with the current-dependent average, E— (°) = EoJ 2 /J 2 rit , and current-dependent dispersion, 

TE J 2 1 Jcrit at J > Jcrit- In the subcritical region the effective temperature = T/(l — J/J cr it) diverges as 

J — » J crit — 0. This singularity can be avoided with due account of the swap processes (see trajectory in figure 1 d) 
that arc important in the vicinity of critical current J rs J C rit- However, this problem is out of scope of this paper. 

It should be mentioned that qualitatively the statistical properties of the "soft" mode in AFM particle are analo- 
gous to those of FM nanoparticle in the presence of spin-polarized current [ ]: Boltzmann distribution and critical 
behaviour of the effective temperature in subcritical region, and Gaussian-like distribution in overcritical regime. 
However, even for one AFM mode the current dependencies of average energy and dispersion are different, as will be 
discussed below. 

From the practical point of view more informative is distribution in full energy, E = E + + E_ , of AFM nanoparticle. 
Corresponding distribution function, f(E), calculated from (12) and (18) as a conventional probability, is shown in 
figure 3. 

Analysis of the curves in figure 3 calculated for different current values shows that in the subcritical regime ( J/ J cr = 
0.1,0.5,0.9) the distirbution function is asymmetric with maximum at 

£max = 77^— tanl^ 1 — Ri T. (20) 

Average energy of the nanoparticle, E av = (E) = + T g , consists of the noisy component only and diverges 
as J — > Jcrit- Energy fluctuation, 



AE ee ^((E-E av y) = (t« + Tj->) ^ (l + (21) 

diverges in the same way (see figure 4 a). Thus, the quality factor, Q — E av /AE diminishes down to 1 as J — > J crit — 0. 
This tendency is quite obvious if one takes into account the noisy source of the energy in the system. On the other 
hand, the presence of singularity in AE shows that thermal noise can play an important role in the transition from 
equilibrium to nonequilibrium steady state in the vicinity of critical current. 
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Figure 4. Current dependence of the average energy, E av , and energy fluctuations, AE in a) subcritical (J < J cr ) and b) 
overcritical (J > J cr ) regimes. Insets show the current dependence of the quality factor, Q = / AE. The energy Eo = 10T. 
Dashed line shows the current-dependence of the "soft" mode energy E^ = Eo J 2 /J 2 T . 

In the overcritical regime the distribution function f(E) has a Gaussian-like shape with the maximum close to the 
average energy E av . Both i? av and AE grow with current (see figure 4 b). However, the main contribution into E av 
arises from the deterministic (low entropy) current-induced rotation of AFM vector, while AE originated from the 
noise slightly intensified by the current. Thus, in this case the quality factor Q is finite at J = J cr and increases with 
current almost linearly (inset in figure 4 b). This behaviour crucially differs from that for FM nanoparticle. Namely, 
in supercritical regime quality factor of FM oscillator, Qfm °c yl - Jcr/ J, vanishes in the close vicinity of J cr . This 
opens a way for the potential applications of AFM nanoparticles as active elements of spintronics devices and as a 
possible alternative to FM nanooscillators. However, further development of the model is necessary to account for 
Joule heating, current fluctuations etc. 

In summary, we considered the current-induced dynamics of AFM nanoparticle in the presence of white Gaussian 
noise which originates from the random magnetic fields. We found the stationary energy distribution functions in two 
regimes: subcritical, when the spin-polarized current is too small to reorient AFM vector from the initial equilibrium 
state, and overcritical, when the spin-polarized current keeps steady rotation of AFM vector. Average energy and 
energy fluctuations in the subcritical regime show the critical behaviour as J — > Jmt- This can be used to facilitate 
the current-induced reorientation of AFM vector. In the overcritical regime the quality factor of AFM particle as 
nanooscillator can be increased by adjusting the current value. 



ACKNOWLEDGMENTS 

The authors are grateful to B.I. Lev for fruitful discussions. The paper was partially supported by the grant from 
the Ministry of Education, Science, Youth and Sport of Ukraine and by the Programme of Fundamental Researches 
of the Department of Physics and Astronomy of National Academy of Sciences of Ukraine. 



[1] J. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996) 
[2] L. Berger, Phys. Rev. B 54, 9353 (1996) 

[3] S. I. Kiselev, J. C. Sankey, I. N. Krivorotov, N. C. Emley, R. J. Schoelkopf, R. A. Buhrman, and D. C. Ralph, Nature 425, 
380 (2003) 



9 



[4] A. A. Tulapurkar, Y. Suzuki, A. Fukushima, H. Kubota, H. Maehara, K. Tsunekawa, D. D. Djayaprawira, N. Watanabe, 

and S. Yuasa, Nature 438, 339 (2005) 

[5] G. Gerhart, E. Bankowski, G. A. Melkov, V. S. Tiberkevich, and A. N. Slavin, Phys. Rev. B 76, 024437 (2007) 

[6] A. Yamaguchi, K. Motoi, A. Hirohata, H. Miyajima, Y. Miyashita, and Y. Sanada, Phys. Rev. B 78, 104401 (2008) 

[7] H. V. Gomonay and V. M. Loktev, Phys. Rev. B 81, 144427 (2010) 

[8] J. Linder, Phys. Rev. B 84, 094404 (2011) 

[9] W. H. Rippard, M. R. Pufall, S. Kaka, S. E. Russek, and T. J. Silva, Phys. Rev. Lett. 92, 027201 (2004) 

[10] A. M. Deac, A. Fukushima, H. Kubota, H. Maehara, Y. Suzukia, S. Yuasa, Y. Nagamine, K. Tsunekawa, D. D. 

Djayaprawira, and N. Watanabe, Nature Physics 4, 803 (2008) 

[11] W. F. Brown, Phys. Rev. 130, 1677 (1963) 

[12] D. M. Apalkov and P. B. Visscher, Phys. Rev. B 72, 180405 (2005) 

[13] V. Tiberkevich, A. Slavin, and J.-V. Kim, Appl. Phys. Lett. 91, 192506 (2007) 

[14] J. Swiebodzinski, A. Chudnovskiy, T. Dunn, and A. Kamenev, Phys. Rev. B 82, 144404 (2010) 

[15] T. Dunn and A. Kamenev, Appl. Phys. Lett. 98, 143109 (2011) 

[16] O. Prokopenko, G. Melkov, E. Bankowski, T. Meitzler, V. Tiberkevich, and A. Slavin, Appl. Phys. Lett. 99, 032507 (2011) 

[17] Y. L. Raikher and V. I. Stepanov, J. Exp. Theor. Phys. 107, 435 (2008) 

[18] B. Ouari and Y. P. Kalmykov, Phys. Rev. B 83, 064406 (2011) 

[19] S. K. Mishra, The European Phys. J. B 78, 65 (2010) 

[20] S. K. Mishra and V. Subrahmanyam, Phys. Rev. B 84, 024429 (2011) 

[21] A. S. Borovik-Romanov, E. G. Rudashevskii, E. A. Turov, and V. G. Shavrov, Physics-Uspekhi 27, 642 (1984) 

[22] A. I. Akhiezer, V. G. Bar'yakhtar, and S. V. Peletminskii, Spin Waves, interscience (wiley) ed., North-Holland Series in 

Low Temperature Physics, Vol. 1 (North-Holland, Amsterdam, 1968) 

[23] H. V. Gomonay and V. M. Loktev, ArXiv e-prints(Jul. 2012), arXiv: 1207. 1997 

[24] T. Dunn, A. L. Chudnovskiy, and A. Kamenev, ArXiv e-prints(Oct. 2011), arXiv:1110.3750 

[25] B. I. Lev and A. D. Kiselev, Phys. Rev. E 82, 031101 (2010) 

[26] I. Baryakhtar and B. Ivanov, Sov. J. Low Temp. Phys. 5, 361 (1979) 

[27] S. I. Denisov, T. V. Lyutyy, P. H&nggi, and K. N. Trohidou, Phys. Rev. B 74, 104406 (2006), 

[28] M. Takahashi and M. Tsunoda, J. of Phys. D: Appl. Phys. 35, 2365 (2002), http : //stacks . iop . org/0022-3727/35/i= 
19/a=307 



